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Wind power prediction is a widely used tool for the large-scale integration of intermittent wind- 
powered generators into power systems. Given the cubic relationship between wind speed and wind 
power, accurate forecasting of wind speed is imperative for the estimation of future wind power 
generation output. This paper presents a performance analysis of short-term wind speed prediction 
techniques based on soft computing models (SCMs) formulated on a backpropagation neural network 
(BPNN), a radial basis function neural network (RBFNN), and an adaptive neuro-fuzzy inference system 
(ANFIS). The forecasting performance of the SCMs is augmented by a similar days (SD) method, which 
considers similar historical weather information corresponding to the forecasting day in order to 
determine similar wind speed days for processing. The test results demonstrate that all evaluated SCMs 
incur some level of performance improvement with the addition of SD pre-processing. As an example, 
the SD+ ANFIS model can provide up to 48% improvement in forecasting accuracy when compared to 
the individual ANFIS model alone. 

© 2012 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Wind energy is considered as one of the fastest growing 
energy resources in the world. Globally, wind energy saw an 
average annual growth rate of 31% over the past five years. By 
2020, 12% of the world’s electricity could feasibly come from wind 


power. Wind represents a clean and sustainable source of energy 
and is in abundant supply. However, the increase in wind power 
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penetration requires a number of issues to be addressed including 
the management of wind power generation variability, market 
integration, interconnection standards, power quality, power 
system stability and reliability, etc., and as such, the power grid 
operators and energy traders face a number of challenges due to 
the substantial growth of increased wind power penetration in 
electricity energy systems [1,2]. Since the wind power produced 
by a wind farm depends critically on the stochastic nature of wind 
speed, unexpected variations of a wind power output may 
increase operating costs for the electricity system. 

The relationship between wind speed and wind power is 
highly nonlinear (i.e., basically cubic), hence, any error in the 
wind speed forecast will actually generate a larger error in wind 
power production. When an entire wind farm is considered, the 
relationship becomes more complex as speed and directional 
components of the wind are used by individual turbines to 
achieve an optimal power output of the wind farm [3]. Corre- 
spondingly, the application of short-term wind speed forecasting 
(30-min to 6-h-ahead) is instrumental in the planning of eco- 
nomic load dispatch and load increment/decrement decisions 
making with respect to the management of a significant amount 
of wind power. 

Several methods are reported in the literature for short-term wind 
speed forecasting such as the persistence method, physical modeling 
approaches, time-series techniques, soft computing methods, etc. 
The persistence method, also known as a ‘Naive Predictor’, is 
generally used as a benchmark for comparing other tools for short- 
term wind speed forecasting. This method simply uses the past hour 
wind speed value as the forecast for the next hour. Wind forecasting 
methods are usually first tested against the persistence method in 
order to baseline its performance [3]. Numerical weather prediction 
(NWP) is a physical modeling approach used in forecasting wind that 
utilizes various weather data and operates by solving complex 
mathematical models [4,5]. Also, a number of time-series forecasting 
models have been successfully applied to short-term wind speed 
forecasting. The time-series based model uses historical data to tune 
the model parameters and error minimization occurs when the 
patterns match historical ones. Some examples of time-series models 
are autoregressive moving average (ARMA), autoregressive integrated 
moving average (ARIMA), fractional ARIMA (fARIMA), seasonal ARIMA 
(sARIMA), single exponential technique, double exponential, grey 
predictors [3,6-10]. These time-series models perform well in the 
areas where the data is low frequency in nature such as weekly 
patterns, but can have difficulties when there are rapid variations and 
high-frequency changes of the target signal, e.g., wind speed. 
Furthermore, wind speed is a result of the complex interactions 
between large-scale forcing mechanism, e.g., pressure and tempera- 
ture gradients, rotation of the earth and local characteristics of the 
surface. These attributes make accurate wind speed forecasting more 
difficult and highly model dependent. 

Soft computing methods are well known for their capabilities 
when dealing with non-linear systems and have garnered sig- 
nificant attention in the area of wind forecasting. Soft computing 
is a generalized term that encompasses the fields of neural 
networks (NNs), fuzzy logic, evolutionary computation, machine 
learning and probability reasoning. An advantage of the soft 
computing model (SCM) is its lower data dependency compared 
to statistical models and can also handle data non-linearity more 
effectively [11,12]. Another advantage of SCMs is that they do not 
require a priori knowledge of the wind data model [13]. However, 
drawbacks of SCMs include the inherent implicit or hidden input/ 
output data relationship and the possibility of overly excessive 
computational requirements [10]. Among SCMs, NNs have been 
widely used in various forecasting applications and numerous 
types of NNs including backpropagation NN (BPNN), probabilistic 
NN, radial basis function NN (RBFNN), self-organizing feature 


maps (SOFM), cascade correlation NN, extended Kalman filter (EKF)- 
based NN, support vector machine (SVM), and adaptive resonance 
theory NN have been previously suggested [14-17]. Some evolu- 
tionary optimization techniques, e.g., particle swarm optimization 
(PSO), genetic algorithms (GA) have been used for updating the 
weight of a neuron while training the NN [18]. Also, fuzzy theory has 
been used for various forecasting applications [19,20] and another 
approach for short-term wind speed forecasting based on the fuzzy 
ARTMAP technique was explored in [21]. 

Approaches using a hybrid intelligent system for wind speed 
forecasting are becoming more popular, namely neuro-fuzzy 
methods. Fuzzy logic and NNs are natural complementary tools. 
A wide comparison among ARMA, NNs, and adaptive neuro-fuzzy 
inference system (ANFIS) models are presented in [22] for fore- 
casting wind power in five forecasting horizons (1 h, 3 h, 6 h, 12 h, 
and 24 h). Li et al. [23] proposed the use of a Bayesian combina- 
tion method using NN models and showed improved wind speed 
forecasting accuracy. Cadenas and Rivera [24] used a hybrid 
ARIMA-NN model for wind speed forecasting and their results 
demonstrated that the forecasting performance of this hybrid 
model for a fixed prediction horizon was significantly better 
compared to that of the ARIMA and the NN models working 
separately. Blonbou [25] presented an adaptive very short-term 
wind power prediction scheme using NNs as predictor along with 
adaptive Bayesian learning and Gaussian process approximation. 
The results showed that the NN predictor performs better than 
the persistence model and the Bayesian framework permits the 
prediction of the interval within which the generated power 
should be observed. Recently, Bhaskar and Singh [26] applied an 
adaptive wavelet neural network (AWNN) based wind speed 
forecast model and utilizes multiresolution decomposition of 
wind speed signal. It was found that the AWNN model has better 
forecasting accuracy and a faster training ability when compared 
to a conventional feed forward neural network (FFNN). Given the 
wide variations of wind forecasting algorithms, a standardized 
protocol for evaluating short-term wind power prediction sys- 
tems was introduced in [27,28] where the authors also provided 
guidelines for using a minimum set of error measures. It was 
emphasized that a rigorous use of data is required as well as the 
fact that both training and testing data sets should be clearly 
defined and separated. Giebel et al. [29] presented a literature 
overview on the state-of-the-art in short-term wind power 
prediction. The authors stated that wind power prediction tool 
is not “plug-and-play” as it highly depends on wind farm site, and 
when installed to a new site, a considerable effort is needed for 
tuning the models based on the characteristics of the local wind 
profile. As explained here, several forecasting methods are avail- 
able for wind speed forecasting, however, there is still a need of 
an efficient and robust forecast tool. 

This paper describes short-term wind speed forecasting 
approaches by considering various SCMs, and the combined 
approach of an SCM and the similar days (SD) method for data 
pre-processing. The forecasting is applied to data obtained from 
the North Cape wind farm located in the Prince Edward Island 
(PEI), Canada. This paper reports the results obtained for 1-h- 
ahead (case-I) and 3-h-ahead (case-II) wind speed forecasts for 
the randomly chosen days of multiple seasons of the year 2010. 
The forecasting horizon in both the cases are 24 h and 72 h. In this 
paper, the following forecasting procedures are analyzed: 


(i) forecasting based on averaging a selected number of similar 
wind speed days corresponding to forecast day, i.e., SD 
method. 

(ii) forecasting based on considering BPNN only. 

(iii) forecasting based on considering RBFNN only. 
(iv) forecasting based on considering ANFIS only. 
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(v) forecasting based on an integrated forecasting model by 
utilizing SD combined with the SCM where the output of 
the SD method (SD) is fed into the SCM in order to obtain the 
final estimate of wind speed. In other words, the wind speed 
forecasts process is based on averaging a selected number of 
similar wind speed days and refining the results through the 
application of SCM. Three different models were assessed 
here: SD+BPNN, SD+RBFNN, and SD+ANFIS. 


Hence, there were a total seven different forecasting models: 
four individual models, namely SD, BPNN, RBFNN, and ANFIS, and 
three combined models, namely SD+BPNN, SD+RBFNN, and 
SD+ANFIS that were assessed for wind speed forecasting accu- 
racy. It is emphasized that though the performance of seven 
different forecasting models are explored in the paper, this paper 
introduces SD+ANFIS as the proposed short-term wind speed 
forecasting model. To demonstrate the effectiveness of the pro- 
posed hybrid SD+ANFIS model, the proposed model is compared 
with all other individual and combined models as well as with the 
methods based on the linear prediction (LP) and NWP. Further- 
more, 1-h-ahead forecasting results were also compared with the 
persistence method, which is considered as a benchmark method. 
The test results demonstrate a satisfactory forecasting performance 
of all the hybrid models (SD+BPNN, SD+RBFNN, and SD+ANFIS) 
considered in this paper, however, a large reduction in forecast error 
of an individual ANFIS model by more than 40% reflects the 
effectiveness of the proposed hybrid model based on SD+ANFIS. 
The forecasting performance of the proposed SD+ANFIS is not only 
more adequate and effective than the individual forecasting perfor- 
mance of ANFIS, but also it performs better than the single use of SD 
method, BPNN, RBFNN, LP, NWP as well as the integrated 
approaches using SD+BPNN and SD+RBFNN models. In general, 
the test results demonstrate that the proposed SD+ANFIS based 
prediction model performs effectively in transforming the given 
numerical weather data into wind speed predictions with a sig- 
nificant reduction in forecasting errors. 

The proposed hybrid SD+ANFIS model for wind speed predic- 
tion is characterized by: 


e SD method based model performs a linear mapping between 
inputs and outputs. The ANFIS model is a hybrid of two 
intelligent system models and combines the low-level compu- 
tational power of an NN with the high-level reasoning cap- 
ability of a fuzzy system thus providing smoothness and 
adaptability. Hence, the combination of SD and ANFIS provides 
a robust approach to wind speed prediction. 

e Ability of modeling the interaction of wind speed and other 
weather parameters such as wind direction and temperature 
in the forecast process. 

e Effective forecasting performance for multiple seasons with 
higher accuracy. 

e The SD method requires a large amount of historical data for 
modeling the input parameters. 

e The ANFIS model has a higher computational complexity. In 
other words, the learning complexity of this model is highly 
affected when inputs or membership functions are increased. 


While there are several papers discussing the wind speed 
forecasting topic, the original contribution of this paper is an 
upgrade of the existing techniques. The fundamental and novel 
contribution is the choice of some input factors based on SD 
technique for better accuracy and proposing a new hybrid fore- 
casting technique based on the combination of SD method and 
SCM(s) in order to translate the practically available numerical 
weather data from a wind farm into wind speed signals. 


The rest of this paper is organized as follows. Section 2 
presents wind speed forecasting by SD method. Soft computing 
models are discussed in Section 3 followed by the numerical 
results and discussion in Section 4. Section 5 includes the final 
conclusions of the paper. 


2. Similar days method for wind speed forecasting 


The basic principle of the SD method is to match historical 
with similar properties to that of the forecast day [30]. In the 
context of wind speed similarity, minimum-distance Euclidean 
criterion is adopted in this paper in order to evaluate the similar 
wind speed days corresponding to the forecast day. 

A proper selection of input variables is important in order to have 
the better performance of the prediction model. Even though many 
exogenous variables can have a good relationship with wind speed, 
it is not feasible or efficient to consider all of them in the forecasting 
process. Some of the major factors affecting the wind speed (WS) are 
wind direction (WD), temperature (T), pressure, humidity, etc. This 
paper considers an interaction between WS and other weather 
parameters such as WD and T in the wind speed forecast process. 
Other influential aspects such as humidity, pressure, etc., are not 
considered as inputs in this paper; however, these aspects will be 
explored in our future work. For the selection of similar days, the 
following Euclidean norm equation adopts WS at time t and t—1, 
and wind direction at t. 


IENI = \/W(AWS,)? + W2(AWS;_1)2-+ W3(AWD,), (1) 
where 

AWS, = WS- WS?" (2) 
AWS (1 = WS;-1—WSPS*, (3) 
AWD, = WD,—WD?P™, (4) 


where IIENII is Euclidean norm; WS; and WD; are WS and WD on the 
forecast day, respectively. WS;_1 is WS at (t—1) on the forecast day. 
ws?" and WD?" are WS and WD in the historical days, respec- 
tively. WS?" is the WS at (t—1) on the historical days. AWS;_1 and 
AWD; are the WS and WD deviation between the forecast day and 
the similar days in past, respectively. W; (i= 1—3) is the weighting 
factor and it is determined by using the method of least squares 
[31,32]. Eq. (1) is an important equation for the selection of similar 
wind speed days based on the historical data. Note that SD 
represents similar wind speed days throughout this paper. We can 
observe from (1) that an interaction of WS and WD is taken into 
consideration during the SD selection process. In order to minimize 
the input variables without compromising the accuracy of the SD 
based forecasting method, it is observed that WS and WD at the 
same hour, and WS at the previous hour are appropriate inputs to 
find the SD. Note that WD is a circular variable. Authors attempted 
to incorporate circularity effect in (1) and found comparatively 
worse performance of the SD method with high forecasting error. 
Literature review also suggest that there is no significant improve- 
ment in wind speed prediction when sine and cosine of wind 
direction is considered [33,34]. 

Similar days are based on the same season. In Fig. 1, the slanted 
portion indicates the ‘forecasting period’ and D3 = 135 days is the 
time framework for the selection of SD corresponding to forecast 
day, i.e., the SD method considers the previous 45 days (e.g., as 
shown in Fig. 1, this data goes back till August 17, 2010) from the day 
before a forecast day (September 30, 2010 at 23:00 h), and 45 days 
before (August 17, 2009) and after (November14, 2009) from the 
forecast day in the previous year (September30, 2009 at23:00 h) for 
the selection of similar days. SD method gives hourly calculation. 
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It means SD method uses 45324 = 3240 hourly data points. 
Among these hourly data, the SD method chooses the lowest 100 
values from (1) using lIENII criterion, i.e., 100 low values are sorted 
out of 3240. Note that the lower the value of |IENIl, better the 
selection of SD. Then, these 100 values are averaged in order to get 
the SD forecast. Authors did sensitivity analysis by choosing 
various data samples as 10, 20, 30, 50, 75, 100, 150, and 200, and 
it was found that the data sample with 100 produces the best 
result. In Fig. 1, September 30, 2010 at 23:00h is chosen as a 
forecast day for illustration purposes only. It should be noted that if 
the forecast day is changed, similar days are selected in the same 
manner. 

An illustration of the selection of hours for the SD criterion is 
shown in Fig. 2(i) and (ii) where 1-h-ahead and 3-h-ahead 
forecasting procedures utilizing the SD method are described. 
For instance, in Fig. 2(ii), similar wind speed data is selected by 
calculating IIENI| from (1), i.e., wind speed data of SD at time t+1 
are obtained first and it is assumed to be the forecasted wind 
speed Din: at time t+ 1. Note that (1) gives an hourly calculation 
of IIENII. SDi: is an average of 100 selected hours of SD data, and 
using it, SD wind speed forecast SD++2 is obtained for time t+2. 
SDi is again an average of 100 selected hours of SD data. In 
other words, a different set of 100 SD data is selected for each 
hour based on (1). This process continues until the target hour, i.e, 
once the forecasted SD data SDi+2 is obtained at t+2, this SD 


D 
— >» 


2010 


eee a 


8/17/2009 


Time "t" on Forecast day (e.g., 9/30/2010 at 23:00) 


11/14/2009 


Previous 
Year 


same time "t" on historical day, 1 year ago, i.e., 9/30/2009 at 23:00 
LA: forecasting period; D=45 days 


Fig. 1. Evaluation of SD based on historical information. 


AN 

SD py 
E a hour 
t t+1 


e.g., WS at 23:00 
on 9/30/2010 
(current hour) 


WS at 0:00 
on 10/01/2010 
(target hour) 


output is used in order to select similar days wind speed data at 
t+3, i.e., SDi+3, thus providing 3-h-ahead SD based wind speed 
forecasts. For example, in Fig. 2(ii), for calculating SDs, i.e., WS 
at 00:00 on 10/01/2010, the input parameters are WS;, WS;_1, and 
WD,, where t is 23:00 on 9/30/2010. To get the value of SDr 4.2, i.e., 
WS at 01:00 on 10/01/2010, the input parameters are SDi41, WSs, 
and WD,.,.1, where t is 23:00 on 9/30/2010. Similarly, to get the 
value of SD, 43, ie., WS at 02:00 on 10/01/2010, the input 
parameters are SD;+2, SDt+1, and WD;+2, where t is 23:00 on 
9/30/2010. Note that the SD calculation considers the available 
wind direction data for time t+1 and t+2. In this paper, SD 
forecasts are then considered as inputs to the SCMs in order to 
produce the final forecasts. The described procedure for the SD 
criterion is shown in Fig. 3. 


Select similar wind speed 
days by calculating ||EN|| 
using WS, WS;.,, WD; 


Wind speed data of SD at time t+1 is 
obtained above and assumed to be 
forecasted SD wind speed, i.e., Daa 


Repeat the above step and get the forecasted 
wind speed at time t+h-1 and use this data to 
obtain final SD wind speed at target hour t+h 


SD forecast output will then be 
considered as input to SCM 


where h is hour-ahead forecasting, i.e., h=1 for 1-hour-ahead 
forecasting and h=3 for 3-hour-ahead forecasting 


Fig. 3. Forecasting process based on SD method. 


where WS: wind speed 


t: current hour 
AN 


SD p4 :Selected similar days 


wind speed data 


(i) Selection of similar wind speed data for 1-hour-ahead forecasting. 


N ~ 
SD t+1 SD u2 


“aN 
SD u3 


t t1 


e.g., WS at 23:00 WS at 00:00 
on 9/30/2010 on 10/01/2010 
(current hour) 


WS at 01:00 
on 10/01/2010 


t+2 t+3 


WS at 02:00 
on 10/01/2010 
(target hour) 


(ii) Selection of similar wind speed data for 3-hour-ahead forecasting. 


Fig. 2. Illustration of 1-h-ahead and 3-h-ahead SD based wind speed forecast. 
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3. Soft computing models for wind speed forecasting 


Due to their strong learning, cognitive ability and good 
tolerance to uncertainty, SCMs are suitable for the wind speed 
forecasting application. In this paper, BPNN, RBFNN and ANFIS are 
the considered SCMs for wind speed forecasting, and these SCMs 
are separately combined with the SD method to achieve better 
forecasting performance. 


3.1. Backpropagation neural network 


The BPNN is one of the most popular forecasting tools among 
various SCMs. The network architecture is shown in Fig. 4. The 
network’s weights and biases are the keys to the network’s ability to 
forecast wind speed accurately. In each iteration of the simulation, 
the weights and biases are adjusted to minimize the error function 
[35,36]. The interconnectivity among nodes of an NN allows for the 
ability to learn relationship between input-output pairs. The 
“learned” relationship between the input and target output is 
developed during the training phase. The network then obtains 
the wind forecast output from the test set using the trained 
network. In this paper, the BPNN is learned by using Levenberg- 
Marquardt algorithm. The following energy function is used as a 
stopping criterion during the learning procedure [31,37]: 


e= 5>_(0-0*Y, (5) 


where O and O* are network output and desired output, respectively. 

The three-layered BPNN model, as shown in Fig. 4, considers 
WS, WD and T at time t, and the SD forecasted values (SD; 41, 
en Sa) as inputs. The final forecasts are 3-h-ahead wind 
speed forecasting. Note that for 1-h-ahead forecasting, the archi- 
tecture of the BPNN is the same as that of Fig. 4 except that it does 
not consider Shas and Deis as inputs. 


3.2. Radial basis function neural network 


The RBFNN model has a form similar to BPNN. The structure of 
the RBFNN model is shown in Fig. 5. The hidden layer neurons in 
RBFNN contain the radial basis function, which differs it from the 


Input layer Hidden layer Output layer 


wind speed 
forecasts 


“NN Wi) 
KI WD 
1 NaN 


WY 


Fig. 4. Architecture of the BPNN model for 3-h-ahead forecasting. 


Input layer Hidden layer Output layer 


wind speed 
forecasts 
— 


PIWA X Y 
NRA 
WAN OG 
AEX 


Fig. 5. Architecture of the RBFNN model for 3-h-ahead forecasting. 


BPNN, and this paper considers the Gaussian function as a radial 
basis function. A function is radial basis if its output depends on 
the distance of the input from a given stored vector [38]. In the 
RBFNN model, the input is sent to the hidden layer without being 
weighted and the final layer output is the weighted sum of the 
hidden layer’s outputs without transformation. Input parameters 
considered in the RBFNN model are similar to that of Fig. 4 in both 
1-h-ahead and 3-h-ahead forecasting. The output of the hidden 
neurons H,, is converted by the radial basis function neuron and 
is expressed as [39] 


Hm(X) =f (pllx;—Cill) = exp(xi—C;)? /0?), (6) 


where m is hidden neuron and ø is the parameter, which controls 
the smoothness of the Gaussian function ¢. x; is the input vector 
and C; refers to the center of point input vector x. Each hidden- 
unit output is obtained by closeness of the input vector x; to an 
n-dimensional parameter vector C; associated with the ith hidden 
unit [39]. The output of the hidden layer is multiplied by a weight 
associated with the neuron and passed to the summation, which 
adds up the weighted value, and it is denoted by 


N 

yi= >> (WHn), (7) 
n=1 

where w, denotes the weight associated hidden layer to output 

layer. 


3.3. Adaptive neuro-fuzzy inference system 


The NNs have a low-level computational structure, whereas 
fuzzy logic deals with reasoning on a higher level. However, the 
disadvantage of using fuzzy logic is that fuzzy systems are unable 
to learn and cannot adjust themselves. Hence, a hybrid network 
based on fuzzy system and NN can be a good strategy to improve 
short-term wind speed forecasting accuracy [37]. This hybrid 
system has thus advantages of both neural networks (e.g., learn- 
ing capabilities, optimization capabilities), and fuzzy systems 
(e.g., human-like ‘if-then’ rules, and ease of incorporating expert 
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knowledge available sometimes in linguistic terms) [40-42]. The 
five-layered ANFIS model is as shown in Fig. 6. Note that the SD 
method is combined with ANFIS model making an integrated 
forecasting model as SD+ANFIS. Among different individual and 
integrated SCMs implemented in the paper, this paper focuses on 
the SD+ANFIS technique as the proposed wind speed forecasting 
model. As shown in Fig. 6, each layer of the ANFIS has its own 
task, which is briefly described as: 

Layer 1 is an input membership function layer. The inputs are 
fuzzified here. Three membership functions are used, and A, B, C, 
D use bell-shaped membership functions [41]. 

Layer 2 is a rule layer and the outputs of which represent the 
membership grade of the corresponding inputs. The node gen- 
erates the output by cross multiplying all the incoming signals. 
Each output node represents the firing strength of a rule. If M is 
the number of the membership function and K is the number of 
the input, then M“ gives the total number of rules. Given the large 
number of rules used in this paper, the connection between layer 
1 and layer 2, and also between layer 2 and layer 3 are shown 
partially in Fig. 6. 

Each node of the layer 3 receives inputs from all nodes from 
the layer 2 and calculates the normalized firing strength, which is 
the ratio of the firing strength of a given rule to the sum of firing 
strengths of all rules [42]. 

Layer 4 is the defuzzification layer. Each neuron in this layer is 
connected to the output of the layer 3 and also connects to all 
inputs. The output of the layer 4 is the function of consequent 
parameters of rules created from layer 3 and input variables [37]. 

Layer 5 is the output layer. The single node in this layer is a 
fixed node, which computes the overall output by summing all 
incoming signals and produces wind speed forecasts. 

A higher number of inputs or membership functions can 
produce a large number of fuzzy rules, which can increase the 
learning complexity. To reduce this complexity, SDi41. SDi +2 and 
SDi43 are averaged, i.e., SDs, t),t;, and then considered as one 
of the input parameters along with WS, WD and T at time t 
(see Fig. 6). 


Layer 1 Layer 2 Layer 3 


3.4. Training and testing procedure for SCMs 


The procedure for the proposed algorithm to forecast h-hour- 
ahead wind speed including learning and forecasting methodol- 
ogy of the SCMs is described below: 

Step-1: Get historical WS, WD and T data. 

Step-2: Select h-hour-ahead (where h=1 and 3) and forecast- 
ing horizon (24 or 72 h). 

Step-3: Select SD for the first learning day of SCM from the time 
framework for the selection of SD (see Fig. 1), which considers the 
previous 45 days from the day before a forecast day, and 45 days 
before and after from the forecast day in the previous year, i.e. 
45+3:24 = 3240 hourly sample data. Among these hourly data, the 
SD method chooses the lowest 100 values from Eq. (1) using IIENII 
criterion. Then, 100 selected SD are averaged to get SD. 

Step-4: Training data of SCM. The online training of the SCM 
has been adopted in this paper. Past 30 days data from the day 
before a forecast day are used for training the SCM. Training data 
includes WS, WD, T, and SD where SD is the SD forecast output. 

Step-5: SCM learning algorithm is applied to train the network. 
The SCM is trained for all the days of learning range (30 days) 
using WS, WD, T including SD obtained from Step-3. Note that 
each learning day of the learning range has one learning data set 
and that makes the total data set of 30%24=720. Once the 
training data set reaches the 720th hour of the learning range, 
this ends the cycle-1, i.e., epoch-1. The stopping criterion is an 
early stopping with maximum number of iterations as 1000 
epochs. Table 1 presents the training parameter of SCMs. 

Step-6: The training of the SCM is completed in Step-5. Then, we 
select the SD wind speed data, which is an average of 100 selected 
SD data, corresponding to the forecast day. This average SD data is 
the SD forecast and then fed to the SCM for the final forecast. 

Step-7: Input variables of the applied SD+SCM includes: 
(i) SDevt (for 1-h ahead forecasting) or Ciawee SDr42, Cae 
(for 3-h-ahead forecasting); (ii) actual wind speed data, WSs; 
(iii) actual wind direction data, WD,, and (iv) actual temperature 
data, T;. 


Layer 4 Layer 5 


wind speed 
forecasts 


Fig. 6. Architecture of the ANFIS model for 3-h-ahead forecasting. 
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Table 1 
Model parameter selection for SCMs. 


Model Parameters 

BPNN Learning rate (7) = 0.9, momentum (a) = 0.2, 
number of neurons in hidden layer=9 (for 6 inputs), 
training function=TRAINLM, 
epoch= 1000 


RBFNN Spread of radial basis functions=0.8, 
maximum number of neurons= 13, 
training function=“‘newrb” Matlab routine, 


epoch= 1000 


ANFIS Number of membership functions (1) = 3, 

input/output membership function type=generalized bell-shaped, 
optimization method = HYBRID (combination of least-squares 
method and BP gradient descent method), 


epoch= 1000 


Step-8: This step involves the final forecasting of hourly wind 
speed using the combined SD+SCM. The input variables 
described in Step-7 are used in the trained SCM and we obtain 
the hourly wind speed forecasts. 

With the change in forecast time, the initial conditions of the 
SCM parameters remain unchanged. If the forecast time is 
changed, SCM is retrained to obtain the relationship among WS, 
WD and T around the forecast day. 


4. Numerical results and discussion 


As stated previously, this paper presents a short-term wind 
speed forecasting technique based on averaging wind speed 
forecast of the SD method and refining the results through the 
application of SCMs. The forecasting approaches were tested 
using the real data of North Cape wind farm located in the PEI, 
Canada. Data sets of WS, WD and T are recorded in 10-min 
interval in the North Cape wind farm. The measured data of the 
six 10-min interval over an hour are averaged to obtain hourly 
data used in this paper. Given hourly wind speed forecasts are 
essential to have an effective output of the wind facilities and 
wind power production forecast. The power system operator 
requires hourly wind power production forecast in advance in 
order to carry out planning operations such as unit commitment 
and economic dispatch. 

For a fair comparison, the same data and weather parameters are 
considered in all the forecasting models, i.e., SD, BPNN, RBFNN, 
ANFIS, SD+BPNN, SD+RBFNN, and SD+ANFIS. In a short form, 
SD+SCMs represent SD+BPNN, SD+RBFNN, and SD+ ANFIS. The 
wind speed forecasting has been carried out using two major cases: 


e Case-I: forecasting horizon considering 1-h-ahead, i.e., fore- 
casting wind speed one hour before real-time for two different 
forecasting look-ahead times of 24 h (daily forecasts) and 72 h 
(three-day forecasts) 

e Case-II: forecasting horizon considering 3-h-ahead, i.e., fore- 
casting wind speed three hours before real-time for two 
different forecasting look-ahead times of 24 h (daily forecasts) 
and 72 h (three-day forecasts). 


It is emphasized that cases-I and -II are carried out by 
randomly choosing one day (daily forecasts) and three continuous 
days (three-day forecasts) from each of the seasons of the year 
2010. The actual flow of the wind speed forecasting process using 
SD and SCMs is shown in Fig. 7 where fi, fo, fs, fa, fs, fe, and fy 
represent forecast from the SD method, BPNN, RBFNN, ANFIS, 
SD+BPNN, SD+RBFNN, and SD+ANFIS models, respectively. 
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Ws, i BPNN 
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H ig M 
| SD+BPNN l 
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WD; P A <x— WS, 
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D 
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(proposed) 


Fig. 7. Flow of wind speed forecasting using an integrated approach of similar 
days and soft computing model(s). 


4.1. Accuracy measures 


For all the days under study, three types of accuracy measures 
are computed to evaluate the performance of the forecasting 
models. The criteria to compare the performance of the proposed 
SD and SCMs are: 


1 24 | wst! WS forecast | 
MAPE gay = 34 5 l ws Tue 
i=1 i 


where MAPEgay is the daily forecast error, WS;™° is the actual wind 
speed in hour i, wee is the predicted wind speed for 
that hour. 

Analogous to the daily error, the three-day forecast error is 
computed as 


100%, (8) 


1 72 ws" -WS forecast 
MAPE three-day ~ 79 5 | - wre | 
i=1 i 


100%. (9) 


Additionally, the root mean square error (RMSE) is also 
computed. This parameter is defined as the square root of the 
average of either 24 (daily) or the 72 (three-day) square differ- 
ences between the predicted values and the true ones. 


24 
RMSE gay = A (ws warty, (10) 
i=1 
1 72 g 
RMSE three-day = 55 Sows? ue ws erecasty2. (11) 
t= 


Furthermore, in order to assess the prediction capacity of the 
forecasting models, mean absolute error (MAE), which provides 
an indication of error ranges is calculated. 


M 
MAE = D Įwste —wshrecast ; aay 


i=1 
where the value of M is 24 and 72 for daily and three-day 
forecasting, respectively. 


4.2. 1-h-ahead wind speed forecasting results 


Table 2 presents the case-I results obtained from all the chosen 
forecasting models and the results are compared with the persis- 
tence method. As seen from Table 2, the MAPE values obtained from 
the persistence method are 13.23%, 18.69%, 16.91%, and 6.75% for 
winter, spring, summer and fall, respectively. Note that December 20 
(Monday), May 21 (Friday), July 5 (Monday) and October 10 (Sunday) 
of the year 2010 have been chosen from the seasons-—winter, spring, 
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Table 2 
Comparison of 1-h-ahead forecasting performance of the chosen models for the 
forecasting horizon of 24 h. 


Season Error Persis SD BPNN RBFNN ANFIS fs fe fF 


Winter MAPE 13.23 14.82 13.88 13.67 12.85 8.89 8.67 7.65 
MAE 1.11 1.32 1.21 1.19 1.15 0.77 0.68 0.58 
RMSE 1.38 1.74 136 1.43 1.34 0.98 1.02 0.73 


Spring MAPE 18.69 13.72 15.82 14.11 13.88 10.88 9.23 8.02 
MAE 0.93 0.52 0.62 0.50 0.56 0.45 0.34 0.35 
RMSE 1.21 0.66 0.74 0.71 0.79 0.54 0.59 0.58 


Summer MAPE 16.91 9.52 14.37 13.92 1413 9.45 9.65 7.27 
MAE 1.01 084 1.28 1.13 1.19 0.99 0.97 0.68 
RMSE 1.34 1.16 1.58 1.45 1.38 1.25 1.28 0.97 


Fall MAPE 6.75 8.72 11.75 10.36 11.06 685 7.46 6.17 
MAE 0.92 0.82 1.02 0.98 0.91 0.87 0.89 0.65 
RMSE 1.17 122 14 1.25 1.32 1.2 1.35 0.94 


Persis: persistence method; f;: SD+BPNN; fẹ: SD+RBFNN; f7: proposed 
SD+ANFIS model. 


summer and fall, respectively. The MAPE values obtained from the 
SD, BPNN, RBFNN and ANFIS models in all the seasons show 
inconsistent results, i.e., these models are found to perform some- 
times better than the persistence method in one season and worse in 
others. In winter, only the ANFIS model (12.85%) shows better 
performance than the other models over the persistence method 
(13.23%). However, the persistence method performs poorer in 
comparison to the SD, BPNN, RBFNN, and ANFIS models in both 
the spring and summer seasons. In fall, the MAPE values from all 
these single models are in the range of around 8%-11% and do not do 
as well as the persistence method (6.75%). From these results, it is 
observed that the performances of the considered SCMs and the SD 
method are seasonal sensitive and they produce inconsistent fore- 
casting results. 

In order to cope with the fluctuated performances of the SCMs 
and SD, this paper presents an integrated model that combines 
the SD method with the SCM, i.e., SD+SCM. With the application 
of SD+ SCM, the errors are greatly improved as it can be observed 
in Table 2 that when SD is combined with ANFIS, the MAPE value 
is found to be lower (7.65%) than that of single use of SD (14.82%) 
and ANFIS (12.85%) in winter. So, with an introduction of SD into 
ANFIS model, the SD+ANFIS output resulted into an improve- 
ment in error by 40.46% over the ANFIS model. The MAPE values 
obtained from the SD+ANFIS demonstrate that there is a large 
improvement in error approximately in the range of 40%-48% in 
all the seasons over the direct use of ANFIS model only. An 
improvement in error from the SD+RBFNN and SD+BPNN 
models is found to be in the range of around 27%-36% and 
31%-41% over the single use of RBFNN and BPNN, respectively. 

The histogram as shown in Fig. 8 compares the MAPE values 
for 1-h-ahead forecasting. Note that Fig. 8 excludes the MAPE 
values obtained from BPNN, RBFNN, SD+BPNN and SD+RBFNN 
models as the focus here is to compare the MAPE obtained from 
the proposed SD+ANFIS model with the persistence, SD and 
ANFIS, and to illustrate a significant improvement in error when 
the ANFIS model is combined with the SD method compared to a 
single use of ANFIS model. 

In order to demonstrate the prediction capabilities of the 
considered models, Table 2 also presents MAE and RMSE values 
for all the test cases. It is observed that the MAE and RMSE errors 
of low values are obtained from the SD+SCM models. In all the 
test cases, the SD method when combined with the SCM outper- 
forms the direct use of a single SCM and compare favorably well 
over the persistence method. The results obtained from the 
proposed SD+ANFIS demonstrates the higher accuracy and 
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Fig. 8. 1-h-ahead errors for forecasting horizon of 24 h. 
Table 3 


Comparison of 1-h-ahead forecasting performance of the chosen models for the 
forecasting horizon of 72 h. 


Season Error Persis SD BPNN RBFNN ANFIS f; fe fF 


Winter MAPE 15.21 15.58 13.65 1455 14.80 10.11 10.55 9,87 
MAE 1.09 1.23 1.16 1.21 1.28 092 0.98 0.84 
RMSE 1.52 1.79 164 1.6 1.53 096 1.35 0.91 


Spring MAPE 22.79 25.69 18.54 17.13 15.67 12.13 12.14 10.1 
MAE 1.10 1.19 1.01 0.82 0.80 057 045 048 
RMSE 1.37 1.33 131 0.90 0.97 086 0.71 0.61 


Summer MAPE 11.73 11.41 15.01 14.93 13.86 12.28 11.3 9.11 
MAE 0.8 0.97 1.24 1.02 1.01 085 0.73 0.56 
RMSE 1.00 1.28 152 1.42 1.12 098 083 0.69 


Fall MAPE 21.08 12.05 14.86 11.77 11.89 12.25 10.71 9.78 
MAE 1.05 1.07 129 1.14 1.10 1.11 1.08 0.84 


Persis: persistence method; fs: SD+BPNN; fs: SD+RBFNN; f}: proposed 
SD+ANFIS model. 


adequacy of this forecasting model. The reason of using 24-h 
prediction by repeating 1-h predictions is to deliver 1-h-ahead 
forecast for the next 24h. Power system operator requires 
different length of forecasting horizon for an operation of inter- 
connected power systems. However, increase in forecasting 
horizon would result into increase in forecasting error. 

This paper also reports the performance of the various tech- 
niques by increasing the forecasting horizon up to 72 h, i.e., three- 
day forecasts for 1-h-ahead wind speed forecasting. Table 3 
presents the three-day-forecasts for the period December 11-13 
(winter), May 17-19 (spring), July 24-26 (summer), and October 
27-29 (fall). The MAPE values using the persistence method are 
obtained in the range of 11.73%-22.79% in all the seasons. As 
expected, the performance of the persistence method is found to 
be deteriorated for the three-day forecasts compared to the daily 
forecast in which the MAPE range was 6.75%-18.69% (Table 2). As 
it can be seen in Table 3, the MAPEs for the SD method, SCMs, and 
SD+SCMs are also found to be increased, but at a lower extent 
when compared to the persistence. The sample average of MAPEs 
obtained from the SD+BPNN and SD+RBFNN models for 1-h- 
ahead (daily forecast) are in the range of 6.85%-10.88% and 
746%-9.65% for all the seasons, respectively. These figures 
increase to only 10.11%-12.25% and 10.71%-12.14%, respectively, 
for the three-day forecasts. Again in the case of increase in 
forecasting horizon, the performance of the proposed SD+ANFIS 
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Table 4 
Comparison of 1-h-ahead monthly forecasting performance. 
Model Error Month 
Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
LP MAPE 21.38 40.73 28.92 31.33 36.08 27.46 26.81 24.19 20.65 19.93 17.82 
MAE 1.12 1.20 0.96 1.34 1.27 1.08 1.11 1.17 0.76 0.84 0.68 
RMSE 1.86 2.34 2.13 2.17 2.34 1.73 2.07 2.25 1.52 1.73 1.32 
NWP MAPE 38.52 49.68 36.52 38.73 39.82 32.09 33.5 37.99 29.55 31.18 26.22 
MAE 1.37 1.22 1.08 1.42 1.33 1.25 1.19 1.36 0.98 1.08 0.94 
RMSE 2.24 2.73 2.33 2.45 2.84 1.98 2.36 2.69 1.88 2.09 1.79 
fi MAPE 15.12 31.85 19.95 21.65 25.20 18.92 20.28 15.76 13.38 13.8 13.26 
MAE 0.71 0.92 0.91 0.99 0.82 0.89 0.71 0.92 0.80 0.76 0.88 
RMSE 0.94 1.18 1.20 1.30 0.99 1.14 0.89 1.16 1.00 0.99 1.19 


fž: proposed SD+ANFIS model. 


forecasting model is superior than the persistence method and 
the other models as well. As it can be seen that the MAPEs 
obtained from the SD+ ANFIS for the three-day forecasts increase 
only by small values (9.11%-10.10%) when compared to that of 
daily forecast (6.17%-8.02%). Furthermore, the low values of MAE 
and RMSE are generated by the proposed SD+ANFIS model and 
the values are within 1 m/s except that the RMSE in fall season is 
1.23 m/s. 

Table 4 shows the two-week forecasts for the period of 
February 16-28, 2010 and monthly forecasts for the period of 
March-December 2010 using LP and NWP models, and the results 
are then compared with the proposed SD+ANFIS. As the total 
data set available in this paper is for the period of January 2009 to 
December 2010, the monthly forecasts for January 2010 and the 
forecasts for February 1-15, 2010 are not reported in this paper. 
The reason behind this is the SD technique considers the data of 
previous 45 days from the day before a forecast day, and 45 days 
before and after from the forecast day in the previous year for the 
selection of similar days. The LP model results presented in 
Table 4 are the forecasts using one-year LP model, the detail of 
which is available in [43]. The NWP model, which is known as 
physical model, results presented in Table 4 are the forecasted 
values obtained from Environment Canada (EC) that uses NWP 
models for their wind speed forecasts. In all the simulated cases 
as shown in Table 4, our proposed model based on SD+ANFIS 
shows better performance compared to that of LP and NWP 
models. Furthermore, to have better comparison, monthly MAPE 
values presented in Table 4 are illustrated in Fig. 9 where we can 
observe the larger values of average monthly MAPEs for LP 
(26.86%) and NWP (35.80%) models compared to our proposed 
model (19.02%). In addition, Fig. 9 also shows the comparison of 
LP, NWP and the proposed model by means of MAPE deviation. 


4.3. 3-h-ahead wind speed forecasting results 


In order to evaluate the prediction capability of the proposed 
SD+ANFIS model, this work also assesses 3-h-ahead forecasting 
(case-II results). In order to forecast and validate the perfor- 
mances of the SD+SCMs, the same days as chosen in case-I have 
been considered. 

Table 5 presents the MAPE results obtained from the SD, each 
SCM and combined SD+SCMs. The sample average daily MAPEs 
obtained from the SD method are 20.91%, 22.22%, 14.21%, and 13.12% 
for winter, spring, summer and fall, respectively. These figures 
increase to 22.52%, 26.11%, 18.10%, and 21.32% (Table 6) when the 
forecasting horizon is 72 h, respectively. As it can be observed in 
Table 5, the forecasting performances of BPNN, RBFNN and ANFIS 
models show inconsistency when compared with the SD method. It 
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Fig. 9. 1-h-ahead MAPE monthly comparison. 


Table 5 
Comparison of 3-h-ahead forecasting performance of the chosen models for the 
forecasting horizon of 24 h. 


Season Error SD BPNN RBFNN ANFIS f; fe f3 
Winter MAPE 20.91 16.12 15.91 14.9 10.53 10.91 9.23 
MAE 2.21 1.64 1.62 1.52 0.96 0.99 0.86 
RMSE 2.89 1.84 1.87 1.69 1.36 1.32 1.05 
Spring MAPE 22.22 20.09 19.23 19.86 14.08 12.00 10.77 
MAE 1.05 1.11 0.92 1.02 0.56 0.47 0.46 
RMSE 1.37 1.31 1.24 1.29 0.7 0.81 0.74 
Summer MAPE 14.21 19.62 18.87 20.17 12.35 11.23 10.34 
MAE 1.25 1.71 1.64 1.72 1.24 1.11 1.09 
RMSE 1.57 1.98 1.96 2.09 1.69 168 1.43 
Fall MAPE 13.12 18.29 18.76 17.94 11.52 11.79 10.36 
MAE 1.31 1.96 2.02 1.86 134 130 1.19 
RMSE 1.69 2.29 2.44 2.09 159 148 1.56 


fs: SD+BPNN; fs: SD+RBFNN; f7: proposed SD+ANFIS model. 


is seen that the BPNN, RBFNN and ANFIS models outperform the SD 
method in winter and spring, but fails to do so in the summer and 
fall. However, when SD is combined with the SCM, the forecasting 
performance is found to be enhanced. As it can be seen in Table 5, 
the MAPE value obtained from the proposed SD+ANFIS is lower 
(9.23%, 10.77%, 10.34% and 10.36%) than the single use of ANFIS 
(14.9%, 19.86%, 20.17% and 17.94%) in all the seasons. In other words, 
by combining SD and ANFIS, the SD+ANFIS output resulted in a 
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Table 6 
Comparison of 3-h-ahead forecasting performance of the chosen models for the 
forecasting horizon of 72 h. 


Season Error SD BPNN RBFNN ANFIS fs fe if 


Winter MAPE 22.52 21.95 20.10 19.97 17.12 15.92 14.88 
MAE 2.51 2.36 2.35 2.34 2.13 1.51 1.37 
RMSE 3.30 2.86 2.82 2.83 2.39 1.92 1.80 


Spring MAPE 26.11 25.94 23.22 22.14 1688 13.88 13.31 
MAE 1.86 1.94 1.68 1.78 0.85 0.51 0.62 
RMSE 2.37 2.15 1.82 1.96 1.21 0.78 0.96 


Summer MAPE 18.1 24.00 20.46 20.36 17.23 15.53 13.84 
MAE 1.58 1.98 1.78 1.83 1.49 1.38 1.29 
RMSE 1.91 2.36 2.78 2.13 1.83 1.72 1.52 


Fall MAPE 21.32 22.82 20.26 18.83 17.79 15.26 14.94 
MAE 2.64 2.75 2.61 2.46 2.07 1.92 1.81 
RMSE 3.02 3.45 2.72 2.71 2.48 2.17 2.04 


fs: SD+BPNN; fg: SD+RBFNN; f7: proposed SD+ANFIS model. 
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Fig. 10. 3-h-ahead errors for forecasting horizon of 72 h. 


38.05%, 45.77%, 48.73% and 42.25% improvement in the MAPE values 
for winter, spring, summer and fall, respectively. 

Table 6 presents 3-h-ahead forecasting results when the 
forecasting horizon is extended to 72 h. Note that in all the test 
cases as shown in Table 6, forecasting errors are significantly 
improved by the proposed SD + ANFIS model compared to the rest 
of the models. The performance of the proposed SD + ANFIS model 
in winter is such that its output (14.88%) resulted in a 25.48% 
improvement in the MAPE value over the single use of ANFIS 
(19.97%), and the percentage error improvements in spring, 
summer and fall are 39.88%, 32.02% and 20.65%, respectively. 
The histogram as shown in Fig. 10 for the three-day forecasts in 
all the seasons illustrates the comparison of MAPE value obtained 
from the proposed SD+ANFIS with the SD method and a single 
use of ANFIS model. 

The values of maximum and minimum percentage error 
improvement obtained by using the SD+ANFIS model are 
48.73% in summer (for 3-h-ahead forecasting when the forecast- 
ing horizon is 24 h) and 17.74% in fall (for 1-h-ahead forecasting 
when the forecasting horizon is 72 h). As demonstrated by the 
test results, it can be summarized that a significant improvement 
in forecasting wind speed errors are achieved from the proposed 
SD+ANFIS technique, providing accurate daily and three-day 
forecasts. 

Table 7 shows errors comparison in the case of larger fore- 
casting horizon. The MAPE values start deteriorating when the 


Table 7 
Performance of the proposed model in different forecasting horizons. 


Season Error SD+ANFIS 


1-h- 3-h- 6-h- 12-h- 24-h- 

ahead ahead ahead ahead ahead 

Winter MAPE 7.65 9.45 13.55 18.99 24.10 
MAE 0.58 0.98 1.53 1.74 1.76 

RMSE 0.73 1.23 1.94 2.19 2.07 

Spring MAPE 8.02 14.22 16.84 18.04 22.62 
MAE 0.35 0.87 1.83 1.96 1.85 

RMSE 0.58 1.04 2.14 2.26 2.31 
Summer MAPE 7.27 9.86 15.28 17.09 21.83 
MAE 0.68 0.89 1.76 1.91 1.94 

RMSE 0.85 0.96 1.86 2.13 2.22 

Fall MAPE 6.17 10.39 14.86 17.91 20.8 
MAE 0.65 0.82 1.55 1.86 1.89 

RMSE 0.82 0.93 1.73 1.97 2.40 


forecasting horizon increases. This is due to the propagation of 
the error. In other words, as SDi: increases, accumulated error 
that depends on forecasting period (hourly, daily, weekly, 
monthly, etc.) will also increase in case of predicting wind speed. 

Forecasting short-term wind speed with a higher rate of 
accuracy is always a major concern for the power system 
operator. This paper contributes to alleviate the major problem 
of short-term wind speed forecasting in which the proposed 
prediction model based on ANFIS combined with the SD method 
is applied to forecast 1-h-ahead and 3-h-ahead wind speed. The 
average computation time required by the proposed hybrid 
SD+ANFIS model for 1-h-ahead daily wind speed forecasts is 
around 1 min using MATLAB on a PC with 4 GB of RAM and a 
2.7-GHz-based processor. 


5. Conclusions 


This paper presented soft computing approaches for predicting 
short-term wind speed. The performances of the BPNN, RBFNN 
and ANFIS models were assessed with and without integration of 
the SD method. Wind speed forecasts are carried out in two 
different forms: (i) 1-h-ahead forecast, and (ii) 3-h-ahead forecast, 
both considering the forecasting horizon of 24 and 72h. The 
interaction between wind speed and other weather parameters 
such as wind direction and temperature were taken into con- 
sideration in the wind speed forecast process. The test results 
obtained for 1-h-ahead and 3-h-ahead forecasts confirm that the 
proposed wind speed forecasting algorithm based on ANFIS 
combined with the SD method is capable of transforming the 
historical numerical weather data into wind speed predictions 
with higher accuracy. Possible future work would be to consider 
other weather variables such as humidity, pressure in order to 
generate wind speed and develop a model for the transformation 
of wind speed to wind power. 
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